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MODELING REACTION CONTROL SYSTEM EFFECTS ON MARS ODYSSEY 


JILL L HANNA*, ZACHARY Q. CHAVIS*, AND RICHARD G WILMOTH* 

Abstract. During the Mars 2001 Odyssey aerobraking mission, NASA Langley Research Center 
performed 6 degree of freedom (6-DOF) simulations to determine rotational motion of the spacecraft. 
The main objective of this study was to assess the reaction control system models and their effects on the 
atmospheric flight of Odyssey. Based on these models, a comparison was made between data derived 
from flight measurements to simulated rotational motion of the spacecraft during aerobraking at Mars. 
The differences between the simulation and flight derived Odyssey data were then used to adjust the 
aerodynamic parameters to achieve a better correlation. 
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1. Introduction. Mars Odyssey successfully completed the aerobraking phase of its mission in 
January of 2002, achieving a nearly circular orbit, and after performing minor orbit corrections, project 
engineers began preparation for the science mapping mission. Successful aerobraking depends on 
numerous factors, not the least of which is maintaining a spacecraft attitude that provides optimum drag 
during the pass through the atmosphere. Although the Odyssey spacecraft was designed to be 
aerodynamically stable and was expected to maintain a specific average orientation with respect to the 
flight path, the small aerodynamic forces and moments exerted by the tenuous atmosphere were 
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insufficient to completely dampen any attitude oscillations induced by targeting alignment errors, 
atmospheric winds, residual inertial motions, etc [1], Therefore, an attitude control system consisting of 
reaction wheels and reaction control system (RCS) thrusters along with appropriate control logic was used 
to maintain the spacecraft attitude within acceptable limits. Control during the drag pass was maintained 
by the RCS thrusters which were also used to desaturate the reaction wheels as necessary and to slew the 
spacecraft as required for solar power collecting and communications during the non-aerobraking portion 
of the orbit. 

RCS thrusters firing in a vacuum or in the low-density Mars atmosphere encountered by Odyssey 
produce plumes which expand rapidly and may impinge on other parts of the spacecraft as well as interact 
with the flow around the spacecraft during aerobraking. These plumes may produce undesirable effects 
on the aerodynamic control effectiveness. Impingement effects were discovered during the termination 
phase of Magellan [2], Analyses performed for Mars Global Surveyor (MGS) showed that certain 
thruster firings at the higher densities expected for MGS aerobraking could actually produce thrust 
reversing [3], However, the Odyssey RCS thrusters were considerably smaller than those used on MGS, 
and such effects were expected to be small. Furthermore, the physical arrangement of the Odyssey RCS 
thrusters is such that the principle thrust directions are nearly perpendicular to the flight direction (see 
Figure 1), and any interactions with the flow around the spacecraft would be considerably different than 
for MGS. Therefore, pre-flight analyses were conducted using a simple engineering plume model 
coupled with flow simulations using the Direct Simulation Monte Carlo (DSMC) method to assess the 
magnitude of both impingement and plume -induced flow interaction effects on Odyssey aerodynamics 
[4,5], The results from these simulations were used to augment the aerodynamic database already 
developed to extract atmospheric densities from accelerometer data and for other trajectory simulations 
used in various mission support tasks [1]. These aerodynamic databases both with and without RCS 
thruster effects were then used in 6-DOF simulations to assess the RCS model by comparisons to flight 
data. 

This paper provides a brief description of the RCS plume model and the DSMC simulations used 
for flow field and aerodynamic predictions. The effects of firing various RCS thrusters will be presented 
for Odyssey aerobraking conditions, and the incorporation of the predicted increments in forces and 
moments into 6-DOF simulations will be discussed. Finally, comparisons of the 6-DOF predictions with 
flight data will be provided for selected aerobraking passes. 


2. ACS description. The Attitude Control System (ACS) control frame for Odyssey is 
composed of two star cameras, two sun sensors, two local Inertial Measurement Units (IMU), two sets of 
two-axis gimbals, and four reaction wheels. 1 The IMUs contain the accelerometers from which 
atmospheric density is derived using Eq. (1). 
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This density profile is input to the POST 6-DOF simulation. The reaction wheel desaturation profile is 
also input to 6-DOF POST. The remaining components of the ACS are not considered here. 


2.1 Reaction wheel model. The Mars Odyssey ACS contained four reaction wheels: three axial 
and one skew, offset in the assembly frame so that equal torque is achieved in the x, y, and z spacecraft 
coordinate frame. 1 The reaction wheels were modeled in 6-DOF POST. Telemetry received from the 
Odyssey spacecraft reported reaction wheel speeds. With the knowledge of reaction wheel gains and 
inertias, the reaction wheel speeds were converted to a net reaction wheel torque and summed into the 
environmental moments. 


1 Mars Surveyor Program ’01 - Orbiter ACS Hardware Coordinate Frame Definitions and Transformations. Algorithm: 
LIB-6. Revision: 7. Revision Date: 2/5/01. 
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2.2 Plume aerodynamic model. The Odyssey RCS consists of four monopropellant hydrazine 
thrusters located approximately as shown in Figure 1. 


RCS-3 

RCS-4 



RCS-2 

RCS-1 


Sketch of 
RCS 
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and 

Thrust 

Vectors 


Bottom View 



Front View 


Fig. 1 . Sketch of Odyssey spacecraft showing RCS locations. 


These thrusters have a nominal thrust of 0.88 N, and their axes are canted such that they can be fired in 
various combinations to provide three-axis attitude control (see Table 1). 


Table 1 

RCS firings that control individual axis rotations. 


Axis 

RCS- 1 

RCS-2 

RCS-3 

RCS-4 

+x 

On 

Off 

Off 

On 

-X 

Off 

On 

On 

Off 

+y 

On 

On 

Off 

Off 

-y 

Off 

Off 

On 

On 

+z 

On 

Off 

On 

Off 

-z 

Off 

On 

Off 

On 


Outside the atmosphere, the thruster plumes expand rapidly and can produce impingement primarily on 
the backside of the solar panels. During aerobraking, the plumes may also interact with the flow around 
the spacecraft thereby altering the aerodynamics. Both of these conditions produce incremental forces 
and moments on the spacecraft in addition to those produced by the thrust of the nozzles themselves that 
must be taken into account for accurate attitude modeling. A complete analysis of these RCS 
aerodynamic effects requires modeling of the flow inside the thruster, the flow in the plume, and the 
resulting flow around the spacecraft. Since the latter effects depend on both the aerodynamic (relative 
wind) attitude of the spacecraft and the atmospheric density, the ultimate product of these analyses is an 
RCS aerodynamic database that can provide these incremental forces and moments throughout the 
aerobraking pass, i.e., as a function of attitude and density. 

2.3. Plume flow model. The internal nozzle flow was computed using a Navier-Stokes solverto 
provide the exit plane properties [6], The computations started upstream of the throat using the stagnation 
chamber pressure of 2.034 MPa and temperature of 1 167 K. The geometry upstream of the throat was 
approximated to provide a smooth convergent section; the throat diameter was 0.29 cm and the divergent 
section was a 15° half-angle cone with an exit-to-throat area ratio of 100. The internal flow was assumed 
to be laminar, and the gas was modeled as a perfect gas. 
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The plume flow model used in the present study was devised by Woronowicz and adapted for 
Odyssey by Chavis [6,4], The model is used to compute properties in the plume based on source-flow 
principles and derived from a free-molecular formulation of conservation of mass, momentum, and 
energy. The derivation assumes that the flow expands radially from each of a distribution of point 
sources whose properties are specified based on the exit plane solution described above. Although the 
free-molecular description of the flow is not valid in the continuum core of the plume and the source flow 
assumption breaks down near the nozzle exit, it has been found that the radial expansion assumption gives 
a reasonably accurate approximation of the spatial variations in plume flow properties at sufficiently large 
distances from the exit. Furthermore, this model has been shown to capture much of the functional 
dependence of these properties on nozzle exit conditions and to produce far-field plume properties that are 
comparable to those predicted by full Navier-Stokes computations [6,4], 

3. DSMC flow simulations. The impingement of the plume on nearby spacecraft surfaces and 
the flow around the spacecraft was modeled using the Direct Simulation Monte Carlo (DSMC) method 
[5], The DSMC computations were performed using the DSMC Analysis Code (DAC) of LeBeau which 
was also used separately to compute the basic aerodynamics of Odyssey [7,1]. DSMC models the 
molecular nature of low-density flows by tracking the motion and collisions of millions of individual 
molecules. These flows typically are not in thermodynamic equilibrium and cannot be modeled with 
conventional continuum methods. Macroscopic properties within the flow and forces and moments on 
surfaces within the flow domain are determined by statistical sampling. DAC can treat complex three- 
dimensional geometries and allows non-uniform inflow boundaries to be defined to simulate plumes or 
other inflows whose properties must be defined by other means. 

The plume effects were introduced into the DSMC computations by defining an iso-density 
boundary from the plume source-flow solution to approximate the boundary between continuum and 
transitional flow. Although this boundary is more commonly chosen based on breakdown criteria 
established by Bird, these criteria are typically based on local density gradients which are not accurately 
captured with simple source-flow solutions [5], Therefore, the iso-density condition was chosen for these 
studies based on a ratio of the local-to-freestream momentum flux of 100, which allows the DSMC 
solution to capture the bulk of the interaction of the plume and freestream flows. A typical iso-density 
boundary surface for Odyssey is shown in Figure 2. Since each of the Odyssey RCS thrusters have 
identical nozzle and plume characteristics, this surface is constructed such that it can be translated and 
rotated to represent any of the four thrusters. For multiple thruster firings, multiple surfaces are used. 
The location and size of these surfaces is such that these surfaces do not overlap, so that any plume-plume 
interactions are fully captured by the DSMC computations. 

4. Flow simulation and aerodynamic results. Figure 3 shows surface contours from a typical 
DSMC solution which illustrate both the directed flux onto the spacecraft from the expanding plume and 
the resulting pressure contours on the surface for a typical aerobraking condition. The thrusters 
designated RCS-2 and RCS-3 produced the largest direct impingement effects on the solar panel even 
though they were further away than RCS-1 and RCS-4 because the thrust axes of RCS-2 and RCS-3 were 
canted toward the panel. Flowever, these impingement effects are not as clearly evident in the pressure 
contours of Figure 3, because the atmospheric flow around the spacecraft during aerobraking interacts 
with the plume causing the impingement effects to be diffused. The plume also tends to act as a shield 
which blocks the atmospheric molecules from reaching the panel surface, so that the net effect on the 
aerodynamics is quite complex. 

DSMC simulations were performed with and without RCS firings for various combinations of 
thrusters, various densities (including zero to capture direct impingement effects in the absence of the 
atmosphere), and various spacecraft attitudes. The aerodynamic forces and moments from these 
simulations were then used to construct an database which provided the increments in forces and 
moments that were needed to perform the 6-DOF simulations. A typical set of increments in moment 
coefficients about the spacecraft mechanical axes is shown in Figure 4 for RCS-1 and RCS-2. These 
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increments do not include the contribution from the actual RCS thrust forces but only those caused by 
impingement and aerodynamic flow interaction. Since these increments are shown here as coefficients 
(C m = M C m /qAL, where q is the dynamic pressure and M cm is the moment about the center of mass), their 
values are singular as the density approaches zero. However, in the implementation of the database, the 
zero density increments are replaced by dimensional quantities that represent the vacuum impingement 
effects only 



RCS-2 Plume Inflow 
Surface Used in DSMC 
Simulations 


Fig 2. Plume inflow model used in DSMC simulations of plume impingement and plume-freestream interactions. 
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Fig 3. Plume impingement flux and pressure contours for RCS-2 firing during aerobraking. Atmospheric density = 100 kg/km 3 . 

Spacecraft attitude at zero pitch and yaw. 
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Fig 4. Increments in aerodynamic moments produced by various thruster firings. Spacecraft attitude at zero pitch and yaw. 

The results shown here represent only a small part of the computations required to construct the database 
used for 6-DOF simulations. Further details of the computations and results are provided in Reference 2. 

5. 6-DOF models. Initial inputs to POST 6-DOF were defined at a given time before periapsis 
passage. Depending on the duration of the drag pass, this time ranged from 600 seconds to 400 seconds 
before the time of periapsis. Inputs to the simulation were time, spacecraft position, velocity, orientation, 
angular rates, and reaction wheel speeds. Profiles of density and reaction wheel speeds throughout the 
drag pass were also included in the simulation. 

Orientation, angular rates, thruster on-times, and reaction wheel speeds were obtained through 
spacecraft telemetry. Spacecraft position and velocity were acquired through .spk files generated by the 
Double Precision TRAJectory program (DPTRAJ) at the Jet Propulsion Laboratory. Thruster on-times 
were obtained from Small Forces Files (SFF) generated onboard the spacecraft. The SFF provides 
cumulative effects of delta V over a specific time interval. 2 During Odyssey operations, this time interval 
was set to +600 seconds from periapsis to incorporate any possible RCS thruster firings within the 
atmosphere. 

Density profiles were 7 second running means of the density obtained through accelerometer 
measurements analyzed by the George Washington University accelerometer team at NASA Langley 
Research Center. A 7 second average was used to reduce noise within the density profile while still 
including density spikes or abnormalities within the atmosphere [8], 


2 

“ Software Interface Specification Small Forces File for Mars Polar Lander, Stardust, Genesis and Mars Odyssey 2001 Orbiter 
Version 1.7 February 24, 2001 Custodian : C.Acton Navigation and Ancillary Information Facility. 
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6. 6-DOF results. The 6-DOF POST was used during the aerobraking phase operations at 
NASA Langley Research Center to validate or adjust the aerodynamics subroutine. Comparisons of body 
angular rates and Euler angles were available for every orbit, and any trends in substandard comparisons 
would possibly lead to changes or corrections in the aerodynamics subroutines. Validations of the 6-DOF 
POST were performed on several pre-aerobraking exoatmospheric situations to ensure that proper 
implementation of the non-aerodynamic models in POST were correct. These situations included the 
second Trajectory Correction Maneuver (TCM) on June 20, 2001 and several orbit passes around Mars 
before Odyssey first skimmed the Mars atmosphere. 

7. TCM2 validation. During the cruise phase of the Odyssey mission, there were four 
scheduled TCM bums to correct trajectory dispersions in preparation for insertion into orbit around Mars. 
During the second of these bums designed as an optimization strategy to minimize total required 
propellent, 3 data were received at NASA Langley to validate the RCS and reaction wheel models in the 6- 
DOF simulation. This data required by the 6-DOF simulation of the TCM2 included thruster 
accumulated pulses, quaternions, angular rates, and reaction wheel speeds. From this data, a profile of 
Euler angles and angular rates were simulated using the 6-DOF POST. 

Figure 5 shows 400 seconds of the TCM during which there were few or no thruster firings. If 
there were thrusters fired in this time period, they were fired for shorter than 0.1 second and did not 
record on the telemetry report. This comparison was made to validate the implementation of the reaction 
wheel model. There were slight variations in the reaction wheel speeds that corresponded to slight 
variations in the angular rates of Odyssey. Most of the high frequency oscillations in the flight data is 
attributed to IMU noise. The rates center about zero with no peaks greater than 0.025 deg/s. The 6-DOF 
POST captured the slight changes in rates and therefore validated this aspect of the simulation. 



Time, sec 

FlG 5. TCM2 validation, no RCS firing: Odyssey body angular rates. 


Figure 6 shows a separate phase of the TCM2 during which there were several firings of the RCS 
thrusters. These firings are from RCS-1 and RCS-2 and are most prevalent in the roll rate, resembling 
steep almost instantaneous decreases in rate. The 6-DOF POST captures these firings at the correct time, 
but there is a small overshoot in the thrust that accumulates with time. The thrusters are simply not 
producing as much thrust as the simulation suggests, posing the possibility of thruster impingement 
effects upon the solar arrays that had not been modeled to date. Although the possibility of thruster 
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Mars Surveyor 2001 Mission Plan Revision B, August 2000 JPL 
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impingement had been examined at the time of this validation, the RCS impingement model was not used 
in this validation. Because of the thrust overshoot after approximately 20 seconds, the pitch rate and yaw 
rate also deviate slightly from the flight data. Nevertheless, this test validates the correct implementation 
of the thruster model without impingement effects. 




Fig 6. 


0 20 40 60 80 100 

Time, sec 

TCM2 validation with RCS firing: Odyssey body angular rates. 


8. Exoatmospheric orbit validation. The first orbital passes of Odyssey around Mars 
provided for exoatmospheric orbital data analysis to again validate models that did not depend on 
atmospheric density. Periapsis altitudes prior to P008 were higher than 150 km above the reference 
ellipsoid of Mars, so the 6-DOF POST was run assuming flight in a vacuum. 
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Figure 7 shows Odyssey angular rates 150 seconds prior to and after periapsis 4. The roll rate 
captures all of the thruster firings in this time interval. The vast majority of all thruster firings during 
Odyssey aerobraking were from RCS-1 and RCS-2 to control the roll rate of the s/c (see Table 1). The 
angular roll rate comparison is good to within 0.001 deg/s. The yaw rate also compares extremely well, 
simulating the flight data to within 0.05 deg/s. The pitch rate in the POST resembles the same pattern as 
the flight data with the exception of one point at approximately 65 seconds before periapsis where there is 
a sharp decrease in pitch rate that the simulation does not capture. This is either a thruster firing that the 
small forces file did not report or more probably, a coupling among the axes that the simulation did not 
model. Note that there is a decrease in the yaw rate at the same time which is simulated in the yaw rate 
but does not couple in the pitch rate simulation. If the simulation had followed the flight data at this 
point, the pitch rate would have matched better throughout this test. 

Figure 8 shows the spacecraft roll, pitch, and yaw from the same periapsis. As expected, since 
the simulation of body angular rates and the flight data are so comparable, the Euler angles are similarly 
close; the roll and yaw angles are within two degrees of the flight data after 300 seconds of simulation, 
and the pitch angle is within eight degrees after the same time period. Note that this larger discrepancy is 
again due to the missed decrease in pitch rate at approximately 65 seconds before periapsis. 



-150 -100 -50 0 50 100 15C 



9. Aerobraking results. Once Odyssey entered the Mars atmosphere, the RCS impingement 
model was implemented into 6-DOF POST. Using this model, comparisons of angular rates and roll, 
pitch, and yaw were generated for every orbit during Odyssey operations. The effect of the RCS 
impingement model is shown in figures 9 and 10 of periapsis 161. The blue dotted line indicates the 
Odyssey flight data. The red line indicates the 6-DOF POST without the RCS impingement model. The 
magenta dashed line indicates the 6-DOF POST with the implementation of the RCS impingement model. 
The maximum density of the 7 second running mean of density for this orbit was 49.8 kg/km 3 , a typical 
maximum density approached by Odyssey. The discontinuities in the red line indicate RCS thruster 
firings simulated in 6-DOF POST without the thruster impingement model. As in every orbit during 
aerobraking, the thruster firings appear mostly in the roll rate. The accumulated thrust overshoot in the 
simulation reaches one degree per second after 100 seconds. When the thruster impingement model is 
added to the simulation, this error decreases by 50% within the same time frame. To correct this error 
further, a multiplier was incorporated into the roll moment coefficient with in the impingement model 
only. Note that although RCS-2 and RCS-3 produced the most impingement effects on the spacecraft, 
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RCS-1 and RCS-2 were fired more frequently to control the roll rate, and therefore only a roll moment 
multiplier was necessary. 

The multiplier for orbit 161 increased the impingement roll moment by a factor of 2.35. By 
decreasing the thrust effect, the roll rate was brought closer to the actual flight data. In figure 9 this 
adjusted impingement is indicated by the green ‘+’ line that is comparable to the blue line indicating the 
flight data. As noted in the lower plots of figure 9 this adjustment in the impingement model did not 
much affect the pitch rate or yaw rate. The simulated yaw rate shows good agreement to the flight data. 
Although the pitch rate simulation does not correlate as well with the flight data, the frequency of the 
pitch rate oscillations is similar to flight. 



■160 -100 - 60 0 to WO 150 





Fig 10. Periapsis 161 euler angles. 
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Figure 10 shows the roll, pitch, and yaw from the same orbit 161. Again, the red line of the 
POST simulation without the impingement model does not correlate well with the flight data. After the 
impingement model is added, this error decreases, but it is not until the inclusion of the 2.35 multiplier 
that the simulation compares to within a few degrees. 

The inclusion of a multiplier on the roll moment coefficient of the RCS impingement model was 
introduced during Odyssey aerobraking operations on orbit 46 to adjust the simulation to better correlate 
to flight data. In time, a trend developed in the magnitude of this multiplier. Figure 1 1 shows the 
variation of roll moment multiplier with maximum density of each orbit since periapsis 46. As shown, 
the multiplier decreases as maximum density increases. 



Max Density, kg/krrP3 

Fig 1 1 . Variation of roll moment multiplier with maximum density after periapsis 46. 

During Odyssey aerobraking, as a general rule, the larger density profiles required more thruster 
firings to stay within the control box. Figure 12 shows the on times for each RCS thruster with 
maximum density. RCS thrusters 1 and 2 fired most often, producing a negative roll about the 
spacecraft. As the spacecraft flew through denser regimes of the Mars atmosphere, the RCS-1 and RCS- 
2 fired more to keep the roll angle under control. From figures 1 1 and 1 2, it can be determined that with 
the increased number of thruster firings, the larger accumulated error in roll rate occurred, requiring a 
larger multiplier on the roll moment coefficient in the RCS impingement model to further reduce the 
effect of the simulated thrust. 




Fig 12. Variation of RCS on-times with maximum density after periapsis 46. 
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10. Conclusion. The 2001 Mars Odyssey successfully completed its aerobraking phase with the 
use of RCS thrusters to maintain attitude control within the Mars atmosphere. Several attitude control 
models including reaction wheel and thruster plume models were incorporated into a 6-DOF simulation to 
accurately simulate the rotational motion of the spacecraft during flight. Upon inspection of the thruster 
orientation and plume effects, impingement on the solar array and interaction with atmospheric molecules 
reduced the effectivity of the RCS thrusters. The 6-DOF analysis supported the research that was 
performed pre-flight, and with the use of adjustments to the RCS impingement model, an accurate 
simulation of Odyssey atmospheric flight was performed. 
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